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O ' We study the rate of true vacuum bubble nucleation numerically for a (j) 

field system coupled to a source of thermal noise. We compare in detail the 
cases of additive and multiplicative noise. We pay special attention to the 
choice of initial field configuration, showing the advantages of a version of the 
quenching technique. We advocate a new method of extracting the nucleation 
^ • time scale that employs the full distribution of nucleation times. Large data 

samples are needed to study the initial state configuration choice and to ex- 
tract nucleation times to good precision. The 1+1 dimensional models afford 
large statistics samples in reasonable running times. We find that for both 
additive and multiplicative models, nucleation time distributions are well fit 
by a waiting time, or gamma, distribution for all parameters studied. The 
nucleation rates are a factor three or more slower for the multiplicative com- 
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pared to the additive models with the same dimensionless parameter choices. 
Both cases lead to high confidence level linear fits of In r vs. T-^ plots, m 
agreement with semiclassical nucleation rate predictions. 

I. INTRODUCTION 

In 1969 Langer gave a first-principles formalism for calculating the decay rates of 
metastable states which was applicable to many types of nucleation and growth processes. 
Several limiting cases were discussed, including those of a small energy difference between 
the two phases and overdamped and underdamped systems. Vacuum stability in scalar 
field theory was discussed by Lee and Wick [0in 1974. In contemporaneous work, the 
cosmological consequences of the decay of metastable vacuum states were investigated by 
Zel'dovich, Kobzarev, and Okun and was further developed by Voloshin, Kobzarev, and 
Okun . These authors investigated the idea of bubbles of a true vacuum phase nucleating 
and expanding in a phase of the false vacuum. In 1981 Guth ^ proposed that sufficient 
supercooling of the false vacuum before its decay could lead to enough exponential growth 
of the scale factor (inflation) to solve the flatness and homogeneity problems. This idea has 
itself undergone significant evolution ever since. 

A formalism for calculating nucleation rates, first using semiclassical methods and then 
including quantum corrections, was worked out by Coleman and Coleman and Callan 
in 1977. It was in these latter works that a Lorentz-invariant account of tunneling was 
given. This technique uses a finite-action classical field configuration in Euclidean space, 
called the bounce. The bounce solution is the solution which minimizes the action. They 
also demonstrated how to treat zero modes in the decay rate calculation as well as how to 
determine the critical size of fluctuations which lead to bubbles that grow. In 1983 Linde 
generalized the study to finite temperatures. 

Calculation of decay rates involve solving the equation of motion of the scalar field. For 
the types of potentials which give rise to metastable vacuum states, the equations of motion 
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are non-linear in the field and, unless severe approximations are made, analytic solutions 
cannot be found. The detailed out-of equilibrium, time dependent aspects of metastable 
vacuum decay can only be studied through numerical techniques. False vacuum decay at 
finite temperature has been studied numerically through the use of Langevin equations 
TI| [|ri|] [|T2|. In these works the stochastic term in the Langevin equation was in the form 



of additive noise, a field independent, white noise driving term. Theoretical work |^ |T4 



|T5| [|16| , however, suggest that the noise term may be more complicated and in general may 



be multiplicative and colored. We apply our study of Langevin dynamics to the range of 
parameters appropriate to the semiclassical regime, where there is reason to believe it makes 



sense 



T0[| ||17|| . This way we can be sure about the conclusions drawn from our new nucleation 
time method and about the comparison between the the additive and multiplicative noise 
models. 

We study numerically the decay of the false vacuum in 1 + 1 dimensions, focusing on 
the question of bubble nucleation rate. We consider in detail the effects of two different 
types of noise terms: additive, where the nucleation rate has been studied before in 
1 + 1 dimensions, and multiplicative, which has not been studied in numerical detail in 
this context. Working in one spacial dimension affords several advantages: the system is 
theoretically simple, with renormalization not an issue, for example |l^; the extraction of 
the nucleation time scale by the new method we advocate requires generation of large data 
samples that can be achieved in reasonable times in a 1 + 1 dimensional system; this simple 
system allows us to explore relatively quickly the influence of initial state configuration choice 
on nucleation time distribution; the good statistics allows us to make a sharp comparison 
between the results generated with our exploratory probe of the multiplicative model and 
those generated with the additive model. 

Given our objective we are immediately confronted by the question of the proper defini- 
tion of nucleation time. In order to make full use of the nucleation time data generated with 
each choice of input parameters (temperature, asymmetry parameter, driving amplitude 
and viscosity coefficient), we propose a new description of nucleation time based on a classic 



waiting time distribution With this description of the distribution of nucleation times 
from the simulation data, we extract characteristic decay times for each parameter set. We 
find that the relationship between nucleation time and temperature, derived from semiclas- 
sical calculations, are reproduced by both Langevin models. We also find that nucleation 
times are systematically a factor of three or more longer in multiplicative noise models than 
in additive models, given the same values of scaled input parameters and the same initial 
state preparation in both models. 

In the next section, we briefly review the ideas of vacuum transition in zero and non- 
zero temperature cases. In Section 3 we describe our numerical methods, including the 
preparation of the initial state of the field. In Section 4 we describe our procedure for fitting 
our data and extracting nucleation rates and present our results. In Section 5 we discuss our 
results and draw conclusions and comment on open questions. Three appendices provide 
further details and supporting data for material presented in the text. 



Consider a self-interacting scalar field theory and suppose the field starts off with a value 
= 0, which corresponds to the metastable minimum of the effective potential (hereafter 
referred to as the potential). Classically, transition to the global minimum at = 00 is 
forbidden, but there is a non-zero probability for quantum tunneling through the barrier 
with bubbles of the true vacuum appearing in the metastable phase. The decay probability 
per unit volume can be expressed by the functional integral [§] 



where Se is the Euclidean action. 



^An exponential decay can be fitted to the tails of our distributions, but this would be misleading, 
as we discuss in Sections 4 and 5 



II. QUICK REVIEW OF NUCLEATION RATES 




(1) 
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S,^j,r,H\if^f^\i^f^vm- (2) 
The equation of motion is 

&^ &\t^ _dV_ 
Ot"^ dx^ d(j) 

Notice the positive sign in front of the term on the right. This is the equation of motion of 
the field in the potential — exemplified in Fig. Now the boundary conditions have 
to be suitably chosen. The field starts off at = in the inverted potential and gets to the 
escape point 0i. In the case of degenerate minima the path integral is dominated by the 



non-trivial minimum energy solution of the equation of motion, the instanton [|T^] [^. For 
non-degenerate minima the non-tivial solution to the equation of motion in Euclidean space 
corresponds to the field starting at the unstable minimum, = 0, at r = — oo, 'bouncing' 
off the escape point at r = and asymptotically approaching = again, as r = +oo, the 
bounce solution. Semi-classically, the path integral in Eqn. |l] gives the tunneling rate per 
unit volume as, 

r = /le-^^, (4) 

where Sb is the action, Eqn. ^ evaluated with the solution of the equation of motion 
with boundary conditions 0(r = — oo,x) = 0(r = +oo,x) = 0. The coefficient A contains 
quantum corrections to the tunneling rate. In most cases it is not possible to solve the 
equation of motion analytically with the appropriate boundary conditions. However there 



are limiting cases [||] |g] where exact analytic solutions can be found. It was shown ||2T[ that 
the spherically symmetric, 0(4), invariant bounce is the one with the least action. Then 
with the boundary conditions for the 0(4) symmetric bounce the general solution for 0(p) 
is of the so-called kink form 

<p = \<po[l-tanh{^^)l (5) 

where p = (r^ + x^)^ and dx is the thickness of the transition region. This solution, valid 
when the difference between the energy density in the true vacuum and that in the false 



vacuum is small 0, is called the thin wall solution and the general form is depicted in Fig. |^. 
As the name of the approximation suggests, there is a sharp transition region (or thin wall) 
between the true and false vacuum phases. The interpretation of Fig. Q is that somewhere 
in space, at p = 0, a bubble of the true vacuum, = (pQ, is formed. Away from the center of 
the bubble there is a transition region, which is very narrow on the thin wall approximation. 
Beyond this transition region, and in particular as p — » oo, the field is in the false vacuum 
state = 0. The parameter R appearing in Eqn. ^ is the radius of a critical size bubble 
which can be easily calculated in the thin wall approximation. 

Tunneling at finite temperature can be studied using all of the above formalism if one 
uses the imaginary time technique of Matsubara ||2^ , exploiting the equivalence between four 



dimensional field theory at finite temperature and Euclidean field theory in four dimensions 
with periodicity or anti-periodicity in one of the dimensions. The period is (3 where (3 = 
1/T. One therefore has the requirement that (f){T,p) = (f){T + /5,p). For sufficiently high 
temperatures this additional condition gives the finite temperature action 

Se = f , (6) 

where 5*3 is the three dimensional Euclidean action and T is the temperature. Eqn. ^ is 
valid at sufficiently high temperatures, higher than the temperature at which the period, (3, 
of the bounce solution is equal the critical radius, R, of a bubble P). The solution to the 
equation of motion with the least action is now 0(3) symmetric and the equation of motion 
can be solved in the thin wall approximation as above. 



III. NUMERICAL STUDIES AND TIME EVOLUTION 

The theory of false vacuum decay outlined in the preceding section made clear that 
nucleation can be studied analytically only in a few limiting cases. Out of equilibrium time 
evolution properties can only be studied in detail using numerical techniques. Nucleation 
has been studied numerically through the use of Langevin equations. Early calculations, for 
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example [|T^], used a Langevin equation for the time evolution of the field based solely on 



phenomenological grounds. Other studies have attempted to derive Langevin type equations 
from quantum field theory in the semiclassical limit. Typically the high frequency modes 
are treated as the thermal bath for the low frequency modes, which is identified as "the 
system" ^ []T3| 



For a scalar field theory in one space and one time dimension, a phenomenological 
Langevin type equation can be written as, 

^-^ + .^--V'W+«^.*.T) (7) 

where ^(x, t) is a random force term originating from the coupling of the field to a thermal 
bath. The thermal noise term, ^ (x, t) , and the viscosity, i], are related through the fluctuation 
dissipation theorem [|1^], 



< ^(x, t)^{x', t') >= 2T7]5{x - x')6{t - t'). (8) 

Since Eqn. ^ is central to our work and is used, but seldom discussed, in the field theory 
literature, we provide details on its implementation and its interpretation in our study 
in Appendix A. The appearance of the viscosity, r], in the equations of motion and its 
relation to the thermal fluctuations through Eqn. |] can be understood by considering the 
example of Brownian motion. A Brownian particle suspended in a fluid is subjected to 
random, microscopic, short-duration forces which, over a longer period of time, change the 
macroscopic velocity of the particle. At the same time the random forces in the direction 
opposite to the velocity of the particle serve to retard its motion. The very forces, therefore, 
that are responsible for the macroscopic motion of the suspended particle also are the ones 
that inhibit the motion, namely give rise to the viscosity of the medium. In order to include 
both effects of the random forces, two independent terms are introduced in the equation of 
motion; they are connected through the fluctuation-dissipation relation. 

In Eqn. ^ the fluctuation term is linear. This is referred to as additive noise. However 
there is no reason why one should not consider a general case where the noise is non-linear. 



For example, a Langevin equation with a non-linear thermal fluctuation term has been 
derived, in a set of approximations |T^, by considering a self interacting quantum scalar 



field. It was shown that, at finite temperature, if the short wavelength modes of the field are 
separated out and the two-loop effective potential is used, then the effective action contains 
imaginary terms which are interpreted as coming from Gaussian integrations over auxiliary 
fields, ^. These auxiliary fields appear in the the equation of motion in the form of a random 
fluctuation term. We adopt a phenomenological version of this result in a 1 -|- 1 dimensional 
multiplicative noise model to discover what features of the nucleation rate may differ from 
those of the additive noise case. 

—Qf2 + ^) = + ^)^(^' T). (9) 

Notice that there is now a non-linear viscosity coefficient. In the high temperature limit the 



fluctuation-dissipation relation, Eqn. is valid |T4 . 

We use the Langevin equations with linear and non-linear noise terms, and with the usual 
fluctuation-dissipation theorem to compare the effects of these types of noise on nucleation. 
This is the first such numerical study, as far as we know. 

We re-scale the equations so that all parameters appearing in them are dimensionless. 
We give the detail of the re-scaling in appendix B. We propagate the solution of the equation 
of motion using a staggered leap-frog method ||26|. We have checked that the results do not 
depend on this method. We also repeated simulation runs on a lattice twice as large as the 
one reported here to check that the results show no finite volume (or in 1-dimension finite 
length) effects. We use a lattice of size L = 100 with a space step-size of 6x = .5 and a time 
step of 6t = .1 and with periodic boundary conditions. The value of the field at the n + 1 
time step in terms of the values at n and n — 1, and at position m on the lattice is given, 
for the additive model, Eqn. |^, by 

0n+l,m = '24>n,m " 4>n-l,m + ^i^(V^0„,m " V^n,m - V'{(j)n,m) + ^m,n), (10) 

where, 



and 

V^0n,m = (0n,m+l - 20„.m + (12) 

Gaussian white noise is used for both the hnear( additive) and non-hnear(multiphcative) 
couphngs to the thermal bath. We relate the noise and viscosity through Equation p. The 
discrete form of this equation is 

2Tri 
oxot 

where the 6x and 6t terms are introduced to compensate for the lack of dimensionality of 
the Kroneker deltas. ^ is given by 



Un - ^ sJt^rn,n (14) 

where Gm,n is a Gaussian random number generator of width one. The delta correlation is 
implemented by using a different random kick at each time and at each point on the lattice. 
By looking at approximately 100 bubble profiles we found that fluctuations which take the 
field at five contiguous lattice sites to the true vacuum value always grew to fill the entire 
lattice. Therefore we chose this criterion to determine when a bubble was formed. 



A. Initial Conditions 

The initial conditions in such a numerical study have to be carefully chosen. We looked 
at nucleation times with random, uncorrelated initial field values distributed on the lattice. 
We also looked at nucleation times with short distance correlated and Gaussian distributed 
initial field values. We found that with the first case there was a long delay time during 
which no bubbles were formed. After this delay time there was a distribution of nucleation 
times rather sharply peaked around some value. With the second set of initial conditions, 
i.e. correlated and Gaussian distributed, we found that this delay time was sharply reduced, 
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the distribution rising quickly and peaking much earher. This is illustrated for samples of 
5,000 bubbles in Fig. |. 

The initial field values, whose correlations and distributions are described by Figs. |^ and 
^ are obtained using a quenching technique 1^^. The effect of this initial state preparation 
on the bubble nucleation times is shown by the data indicated by the histogram on the left 
in Fig. 1^. We scatter random field values on the lattice and propagate the solution of the 
equation of motion in a potential with only one minimum. We ensure that the curvature 
at the bottom of the potential closely matches the curvature of the asymmetric double well 
of interest. After some time the field acquires the thermal distribution with the desired 
short distance correlation. At this point, which is t = in the simulation, the potential 
is quenched, i.e., the asymmetric potential is switched on. Nucleation times are recorded 
with reference to this time. Fig. ^ shows the distribution of initial conditions just before 
quenching for both additive noise and multiplicative noise. 

Summarizing, we find that when the initial state is not prepared by the quenching method 
just described, and random, uncorrelated, initial field values are placed on the lattice, there 
is a long delay time for the first bubble to appear, after which a peaked distribution of times 
develops. The delay time apparently corresponds to the time needed for the field to become 
correlated. Preparing the initial state so that it is correlated as shown in Figs. ^ and ^, 
produces the change in the distributions illustrated in Fig. ^ We use correlated inital states 
prepared by the quenching technique in all of the data shown and discussed in the following 
sections of the paper. 



B. Time Evolution 

Fig. ^ shows a series of frames which are snapshots of the lattice at different times 
with additive noise. All results that follow pertain to solutions of the dynamical equations 
rescaled to dimensionless form, as described in Appenxix B. Each point plotted is separated 
by Ax = 1. In the first two frames we see the field throughout the lattice fluctuating about 
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the false vacuum at = 0, with prehminary indication of the formation of a bubble in the 
region between 60 and 70 in the first two frames. The third frame shows an established, 
growing bubble. The fourth frame shows a completely formed bubble. Note that the small 
upward fluctuations in the first frame in the regions around 5 and 95 are gone in the second 
frame. At the top of the profile the field fluctuates about the true vacuum value of = 5.0. 
The expected "tan/i" profile is discernible.0 The matching of the profile at the boundaries 
is due to periodic boundary conditions. In Fig. ^ we display the growth properties of 
bubbles using an example with multiphcative noise. When a fluctuation grows to occupy 
five contiguous lattice sites, we record that time as the nucleation time. 



IV. RESULTS 

A. Identifying Nucleation Times 

The distribution of nucleation times for 5000 bubbles is shown in Fig. ^ at successively 
lower temperatures. We see that the distribution is not symmetric; it rises quickly and there 
is a long tail, which is especially pronounced at low temperatures. From these distributions 
one has to extract nucleation times. We propose that the asymmetric distribution of times 
should be used as a whole. Namely, we seek an appropriate function to fit the entire 
distribution of data at each temperature and asymmetry parameter choice.^ We observe 



^A numerical study of the rate of expansion of the growing bubble can be found in Ref. |Tz 

^Based on remarks in the literature |1C| [11| one might expect that the distribution of nucleation 
times would be exponential. Our results do not show this behavior. Unfortunately, the details 
of the initial field configuration and the actual distributions of nucleation times are not shown in 
these references. In the only paper we know of that does show a nucleation time distribution, for a 
limited statistics sample in a 2 + 1 dimensional model, the shape is in qualitative agreement with 
ours |29|. These authors report exponential fits to the tails of their distributions. This procedure 
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that the distributions rises to their maxima relatively quickly compared to the time over 
which the tail of the distribution persists. We guess that a power law function of time will 
fit this part of the distribution. At long times, a damping function that "overpowers a power 
law" is the falling exponential. The simplest function that has both of the above properties 
and has, at least naively, an appropriate interpretation in our problem, is the classic "waiting 
time" form 

^6t = Ki'-re-^, (15) 
dt ' ^ ^ 

where 5t is the bin size of the histograms. The parameter r is what we use for the nucleation 
time. Both a and r were obtained from fits to the distributions. The normalization, K, 
can be obtained from the requirement that the integral of Eq. [l^ gives the number of 
bubbles that were used to make the histogram. Fig. p| shows the least squares fit to this 



function using a typical additive noise case and Fig. |T0] shows the fit for a multiplicative 
noise case. Figs. ^ and |12] show two more fits with different parameters. These fits were 
made by setting K = 1. From the normalization described above, we checked that this 
value is approximately correct in all cases. For the same parameters in the models, and at 
the same temperature, the nucleation times with multiplicative noise is always longer than 
that for the additive noise case. We emphasize that casting the equations of motion into 
dimensionless form, detailed in Appendix B, makes direct comparison between the models 
possible. 



B. Time vs l/T 

At finite temperature the nucleation rate per unit length, in the semi-classical approxi- 
mation, is given by 



does not give the correct value for the nucleation time because the region of the distribution that 
is fit is not at long enough times to be in the true exponential tail. We develop this point further 
in the work below. 
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r = Ae^ (16) 

where 5'^; is the Euchdean action corresponding to the bounce solution of the equation of 
motion and T is the temperature. A simple manipulation of this equation gives 

lnr = ln{TL)-^ = -ln{AL) + (17) 

where L is the length of the lattice. We extracted the nucleation times from the fits as 
described above and looked at the effects of temperature on nucleation times. To study the 
fluctuations in the parameters in the fits we generated several data sets for a single set of 
parameters. We found little variation in the parameters among runs. The tables in Appendix 
C show a set of fit parameters from different runs. The fluctuations in the parameters a and 
T are small, being less that two percent of the central value. The value of a is especially 
stable from one run to another. We emphasize that, when normalized to 1, the average of 
the fit function is r[a+ 1] x r. This should agree with the average of the data; within (small) 
errors, it does in every case. 

We also looked at the fluctuations in the fits at a lower temperature and found that all 
the fluctuations in the parameters were within five percent of the central value. In fact most 
were within two percent, with a showing smaller variation than r, as in the examples shown 
in the tables in Appendix C. 

Using a global error estimate of five percent in the parameter r to cover all of the cases we 
studied, we fitted straight lines to the data plots of nucleation times vs inverse temperature. 
We show the results for the additive model in Fig. |T3|. The confidence levels of the fits 



are 85%, 98% and 99% for a = .70, .74 and .80 respectively. In Fig. |1^. we show the 
corresponding result for the multiplicative noise model. The confidence levels of the fits are 
95%, 98% and 99% for a = .70, .74 and .80 respectively. Both plots show strong evidence 
that the relationship between Inr and 1/T is linear, reproducing the semiclassical description 
of the finite temperature quantum system, Eqn. 0. 

We see that for both additive and multiplicative noise, the larger values of the asymmetry 
parameter drive nucleation at a faster rate. This is in agreement with nucleation theory: 
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the action for the bounce, which appears in the exponential of the decay rate, is inversely 
proportional to the energy difference between the true and false vacua. The plots also show 
that for a given asymmetry parameter the nucleation time is shorter at higher temperatures 
as one would expect. As mentioned above, the plots also indicate the linear relationship 
between In r and 1/T. The new twist is that t is extracted by our analysis from the full 
distribution as a "waiting time". For long enough times, the waiting time distribution, 
Eq. 16, has an exponential tail with decay constant r. But this is in the time regime 
t^ Tax Init/r), far beyond the point where we run out of simulated events. It would take 
enormous statistics to probe this regime. 



V. DISCUSSION AND CONCLUSION 

We have studied false vacuum decay in 1 + 1 dimensions using Langevin equations on 
a lattice. We looked at cases where the stochastic term in the equation was in the form of 



additive noise, previously studied in [|T0| in 1 + 1, and in the form of the theoretically mo- 
tivated multiplicative noise case |T^. The cosmological applications of multiplicative noise 
have been discussed in ^8|. In past work |[10| [Til decay times were fit by exponen- 



tial distributions. We found distributions of nucleation times that were not exponential, as 
described above. We found that a F, or waiting times, fit with effective multiple incidence 
of about 6, gives a good description of the whole distribution of nucleation times for every 
case studied. In principle, the tail of the fitting function, Eq. 16, is an exponential with 
decay time r. In practice, the time must satisfy t ^ ra x ln{t/T), which far beyond the end 
of our statistics, even with 5, 000 bubble samples. At shorter times, well past the peak of 
the distribution but where data is still sufficient in the falling tail, we can find reasonable, 
approximate exponential fits to the data. However, these fits to the exponential tails do 
not give the same nucleation times as the F distribution, in fact the tails underestimate 
the nucleation times. Fitting the tails of the distribution also requires one to discard most 
of the data and hence potential valuable information on the dynamics of nucleation. We 
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found that for both the additive noise and multiphcative noise Langevin models, the semi- 
classical relationship between nucleation time and temperature is reproduced in the range of 
parameters studied. For the same set of parameters the multiplicative noise equations give 
nucleation times that are longer than in the additive case by factors of three or more. To see 
why this might be so, note that the location of the barrier of the potential, and consequently 
its stable minimum, are at scaled field values significantly greater than one. So apparently 
the quadratic field dependence of the viscosity term in the multiplicative case significantly 
inhibits the thermal hopping process, more than compensating for the enhancement from 
the linear field dependence of the driving term. We tested this idea by varying the value the 
field had to attain for a bubble to be declared formed. We expect that the fractional increase 
in nucleation time as the field value increases should be greater in the multiplicative noise 
case than in the additive. We looked at three different field value requirements, = 5.0, 4.0, 
and 3.0. Tables | and H show that for the same increase in the field value criterion, the 
nucleation time increases by about 4% in the additive noise model and 12% in the multi- 
plicative noise model. While this is not a proof of our conjecture, the trend demonstrates 
the importance of the field dependent viscosity in slowing down the nucleation rate. 

The r distribution which fits the data is the waiting time distribution for a fixed 
number(> 1) of events that obey Poisson statistics. The fit parameter a ~ 6 is fairly 
constant throughout the range of input parameters studied, and this points to some quasi- 
universal feature in the models. The observation that it is of the order of the number of 
lattice sites, 5, that determines our bubble formation criterion, which is consistent with 
the field correlation of the initial configuration, may be a clue to the interpretation of the 
insensitivity of the value of a to the choice of model and of the input parameters for a 
given model. The full interpretation of the successful application of the F distribution to 
this problem and the significance of the approximate universality of the parameter a are 
currently under study. 



VI. APPENDIX A 
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THE FLUCTUATION DISSIPATION THEOREM 



The fluctuation-dissipation theorem ||2|] |25], as the name imphes, relates the stochastic 



fluctuations of a system to its dissipative or irreversible properties. As described earlier its 
application ranges from classical phenomena like Brownian motion to the far more compli- 
cated case of a scalar field coupled to a thermal bath. For a general discussion see ^ 



In numerical studies like the one reported here, the Langevin equation on its own does not 
completely describe the coupling of the field to the thermal bath. The fluctuation-dissipation 
relation, Eqn. ||, is used to connect the stochastic and dissipation terms and hence ensure a 
realistic simulation if the thermal system. The delta function correlations are implemented 
numerically by sampling the distribution of fluctuations at each lattice space site and at 
each time step in propagating the solution of the equation of motion. The exact form of the 
noise distribution and, in particular, its dependence on the viscosity parameter is impor- 
tant because from this we obtain the correct numerical value of temperature. We describe 
here in some detail our implementation of the fluctuation dissipation theorem and the noise 
distribution we sampled in the study. 

The probability distribution of the noise is assumed to be Gaussian of the form, 

dP[Q= ^\ D^e-^^J (18) 



For convenience define 6 = \f2/r\T and = 1/\/2ti6. Then, 

< e(x,t)e(x',t') >=N j Diax,mx\t')e~'^^'-"''"^. (19) 
Introduce a source for the noise, J{x,t), in the exponent. 

< e(x, mx', t')> = N j i{x, mx', f) X 

J dx"dt"^i^!-^p^ + J dx"dt"J{x'\t")i{x",t") ^20) 

J ^ J dx"dt"j{x",t")^{x",t") ^22) 
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Now re-define to be x, t) — ^{x, t) — 9'^J{x, t). Then, 



-Udx"dt"i^ 



Therefore, 



< «x, m.', f)>=Nj Di^^^ X (24) 

g-i / + f / dx"dt"j(x",t")2 ^25) 

= 7V^/2^^ ,,/ (26) 

5J{x,t)6J{x',t') ^ ' 

^ NV2^e'—^J(x',t')e'^I'^"''"'^^"''"^\ (27) 
In the hmit J{x, t) ^ we get, after taking the next functional derivative, 

< ^{x, t)^{x', t') >= NV2^e^5{x - x')S{t - t') = 2r]TS{x - x')S{t - t') (28) 



VII. APPENDIX B 
RE-SCALING THE EQUATIONS OF MOTION 

The parameters rj and ^ have different units in the additive and multiphcative noise 
equations of motion. To compare the results of our numerical study we rescale the equations 
of motion so that all quantities appearing in them are dimensionless. 

Consider the potential, 

V{<P) = - f 0^ + ^<P'. (29) 
V{(f)) appears in the Hamiltonian density, H. 

n^ln' + liVd^r + Vid^) (30) 
The Hamiltionian, H, which has the dimension of energy is obtained from 7i as, 

H = JdxH. (31) 
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V{(f)) therefore has dimension of EL ^. Prom the term (V0)^ in H we get [0] = {EL)^ 
where [•] means "dimension of" . We also have the following; 

[m] = L-^ 
[A] = E-^L-^ 
[a\ ^ E 2L 2^ 

k] = L-\ 

(32) 

The subscript a refers to additive noise. For the multiplicative noise case 77^ and have 
the dimensions, 

[Vm] = E-'L-\ 

[U] = L-\ (33) 

We now rescale the field and parameters appearing in the equation of motion. For the 
purpose of discussion let 0o be the value of the field in the true vacuum.^ When the thin 
wall approximation is not applicable the field value at the escape point at the right of the 
barrier is the more important quantity. For the additive noise case we have, 



(l)Qrn? dt^ (f)oni? dx^ (f)om'^ dt 0o (f)^ rn? 0o 0o"^ 

We now make the following definitions, 

. t 

Va = — , 

m 

k2 



A 



9 : 



^It can also be chosen to be the value that the field takes at the top of the barrier of the potential 
or the value of the field at a zero of the potential. 
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X = m X, 
t = mt, 

T = ^, (35) 

m0Q 

where T follows from the fluctuation-dissipation relation, Eqn. |. The rescaled equation 
of motion is, 

We rescale the multiplicative noise case in a similar manner but the dimensionless vis- 
cosity and noise parameters are, 

= ^ (37) 

m 



and 



L = H (38) 



The dimensionless equation of motion is then 

'dW~d^ ^ 'dt ^ - + ^^0. (39) 

In terms of the original parameters, the pre-scaled T has the dimension of energy but it is 
dimensionless after re-scaling. The key point is that the parameters m, a and A in the scalar 
potential can be used to rescale both equations to dimensionless forms. We have used m 



and 00 = (o; + \/a^^-4m?X)/2 to illustrate the rescaling but other quantities of the same 
dimensions can be used as noted. All quantitive work is done with Eqns. |3^ and ^ so the 
values assigned to a, A determine the value of in the true vacuum in each case. With this 
understanding, the tilde notation is not used in the text when referring to results from the 
rescaled equations. Now that the parameters in both equations of motion are dimensionless, 
we can make direct comparisons of the two models. 
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VIII. APPENDIX C 
EXAMPLE DATA AND FITS USED TO DETERMINE FLUCTUATIONS IN r 

AND a 



Tables |T| through |V11J| and accompanying Figures 15 through 20 show examples of data 
and fits that we used to study sample-to-sample fluctuations. Error estimates used to assign 
confidence levels to linear fits to the In (r) vs. T^^ plots, Figs 13 and 14, are based on these 
studies. 
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FIGURES 
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FIG. 2. The bubble profile, 4>{p), in the thin wall approximation, dx is the narrow transition 
region or wall 



0.25 
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0.05 



-0.05 



FIG. 3. Field correlation < (j){xl)4>{x2) > just before quenching for different temperatures with 
o's representing the highest temperature and A's the lowest. 





(b) 

FIG. 4. Gaussian distribution of initial conditions obtained using the quenching technique, 
^^^^ vs (j), for (a) additive noise and (b) multiplicative noise. The points are data and the curves 
are fits to Gaussians. Here the area under the curve between cp and (p + 6(p is the fraction of lattice 
points having the field value in that range. 
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t 

FIG. 5. Distribution of nucleation times with quenching, left, and with random initial condi- 
tions, right 
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FIG. 6. A series of snapshots of the lattice, (f>{x) vs x, showing the formation of a bubble and 
its subsequent growth for an additive noise system. The absissa is the spatial position and the 
ordinate is the scaled field value. 
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FIG. 7. A series of snapshots of the lattice, 
bubble for a multiplicative noise system. 
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{t = 1500) 

(p(x) vs X, showing the growth of a fully formed 
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FIG. 9. Fit to Eqn. 15 for an additive noise case with T = 1.5 and a = 0.8. From the fit we 



obtained a = 6.8 and r = 26.2 
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FIG. 10. Fit to Eqn. 15 for a multiphcative noise case with T = 1.5 and a = 0.8. From the fit 
we obtained a = 6.7 and r = 84.5 
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FIG. 11. Fit to Eqn. 15 for an additive noise case with T = 1.2 and a = 0.70. 
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FIG. 12. Fit to Eqn. 15 for a multiplicative noise case with T = 1.3 and a = 0.74. 
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FIG. 13. A log-linear plot of r vs. T"^ for additive noise with different asymmetry parameters, 

a=.70, .74 and .80 from the upper to the lower curve. The fits are made with 5% error bars. The 

confidence level of the fits are 85%, 98% and 99% respectively. 
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FIG. 14. A log-linear plot of r vs. for multiplicative noise with different asymmetry 

parameters, a=.70, .74 and .80 from the upper to the lower curve. The fits are made with 5% error 

bars. The confidence level of the fits are 95%, 98% and 99% respectively. 
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FIG. 15. Fit to Eqn. 15 for a multiplicative noise case with T = 1.5 and a = 0.8 
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FIG. 16. Fit to Eqn. 15 for a multiplicative noise case with T = 1.5 and q = 0. 
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250 500 750 1000 1250 1500 1750 
FIG. 17. Fit to Eqn. 15 for a multiplicative noise case with T = 1.5 and a = 0.8. 
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FIG. 19. Fit to Eqn. 15 for a multiplicative noise case with T = 1.5 and a = 0.8. 
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FIG. 20. Fit to Eqn. 15 for a multiplicative noise case with T = 1.5 and a = 0. 
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TABLES 



(j) 5.0 4.0 3.0 

T 27.U2±.18 25.y±.2 24.8±.2 

a 6.63±.005 6.839±.004 6.545±.004 

Tg^g 212±1 204±1 196±1 

TABLE I. Parameters from fits for additive noise with T = 1.5 and a = .80, and A = .1. Each 

column represents a different (j) value which the field was required to reach before the nucleation 

time was recorded. 



(j) 5.U 4.0 3.U 

T 84.24±.43 73.y4±.y 63.8±.6 

a 6.74±.005 6.796±.005 6.796±.004 

Tg^, 654±3 581±3 510±3 

TABLE II. Parameters from fits for multiplicative noise with T = 1.5 and a = .80, and 

A = .1. Each column represents a different (j) value which the field was required to reach before the 

nucleation time was recorded. 
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a 


T 


^avg 


6.83 


26.90 


213 


6.84 


26.71 


211 


6.84 


27.03 


212 


6.83 


26.94 


212 


6.84 


26.99 


213 


O.OO 


Z ( .AA 


ZiO 


D.OO 


OT 1 1 
Z / .ii 


01 Q 

zlo 


D.OO 


Z 1 .60 


zio 


6.83 


26.94 


212 




average 




6.83 


27.02 


212 




cr 




.005 


0.18 


1 



TABLE III. Parameters from fits for additive noise with T = 1.5 and a = .80 
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a 


T 


'^avg 


6.51 


48.27 


366 


6.52 


45.04 


359 


6.51 




362 


6.51 


A A n n 

44.66 


362 


6.51 


A A r* -x 

44.61 


363 


6.51 


A A n A 

44.64 


362 


D.ol 


44. OD 


ODO 


D.Oi 


A A 

44. DO 


ODZ 


o.oi 


A A P.P. 
44.00 


ODZ 


6.51 


44.62 


363 




average 




6.51 


45.03 


362 




(7 




0.003 


1.15 


2 



TABLE IV. Parameters from fits for additive noise with T = 1.5 and a = .74 
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a 


r 




6.04 


103.62 


797 


6.05 


106.59 


787 


6.04 


106.59 


784 


6.05 


103.16 


791 


D.Uo 


iUo.iU 


/ /D 


o.Uo 


iU4. 




D.Uo 


1 no '77 


/Do 


6.04 


104.85 


786 




average 




6.05 


104.67 


784 




a 




0.005 


1.44 


9 



TABLE V. Parameters from fits for additive noise with T = 1.5 and a = .70 
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a 


T 


'^avg 


6.74 


84.48 


f^r A 

654 


6.74 


83.91 


651 


n 'T A 

6.74 


84.58 


657 


6.73 


OA r\o 

84.98 


r* r* r\ 

660 


6.73 


84.23 


654 


6.73 


OA n A 

84.64 


654 


fi 7/1 
D. /4 


00.40 


DOU 


D. /4 


Q/l 1 Q 


D04 


7/1 
0.(4 


Q/i no 


DOZ 


6.74 


83.99 


654 




average 




6.74 


84.24 


654 




(T 




0.005 


0.43 


3 



TABLE VI. Parameters from fits for multiplicative noise with T = 1.5 and a = .80 
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a 


T 


^avg 


6.31 


171.50 


1303 


6.31 


166.91 


1010 

1313 


6.30 


173.00 


1330 


6.32 


168.67 


1300 


6.31 


173.78 


1309 


6.30 


171.67 


1310 


D.oi 


i / O.iD 


ioUo 


D.oi 


1 flQ /I 7 

iDo.4 / 


ioUD 


o.oi 


1 71 1/1 
i / i.i4 


1 Q 1 
ioiO 


6.31 


171.72 


1000 

1323 




average 




6.309 


171.002 


1311 




(7 




0.006 


2.27 


9 



TABLE VII. Parameters from fits for multiplicative noise with T = 1.5 and a = .74 
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a 


r 


^avg 


6.08 


504.00 


3485 


6.06 


509.00 


3754 


6.09 


500.79 


3451 


6.07 


511.79 


3513 


6.09 


495.37 


3495 


6.09 


490.05 


3493 


D.Uo 


A CiP fin 
4yD.DU 


o4o/ 


D.U / 


oUo.oi 


ooiU 


b.U / 


oio. 


Q c;nQ 


6.07 


520.18 


3542 




average 




6.08 


504.70 


3524 




a 




0.011 


9.27 


84 



TABLE VIII. Parameters from fits for multiplicative noise with T = 1.5 and a = .70 
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